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Abstract: In finite-temperature field theory, the cyclic Wilson loop is defined as a rectan- 
gular Wilson loop spanning the whole compactified time direction. In a generic non-abelian 
gauge theory, we calculate the perturbative expansion of the cyclic Wilson loop up to or- 
der g'^. At this order and after charge renormalization, the cyclic Wilson loop is known 
to be ultraviolet divergent. We show that the divergence is not associated with cusps in 
the contour but is instead due to the contour intersecting itself because of the periodic 
boundary conditions. One consequence of this is that the cyclic Wilson loop mixes under 
renormalization with the correlator of two Polyakov loops. The resulting renormalization 
equation is tested up to order and used to resum the leading logarithms associated with 
the intersection divergence. Implications for lattice studies of this operator, which may be 
relevant for the phenomenology of quarkonium at finite temperature, are discussed. 
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1 Introduction 

Wilson loops, ever since their introduction [1], have been an important tool in the study of 
non-abelian gauge theories. In particular, rectangular loops with one spatial and one time 
direction, also known as static Wilson loops, have played a crucial role in the understanding 
of QCD, for their Coulomb-like behaviour at short distances is a manifestation of asymp- 
totic freedom and their area-law behaviour at large distances of confinement. Since static 
Wilson loops in the large time limit are related to the energy of a static quark-antiquark 
pair [2, 3], static Wilson loops are also important quantities for quarkonium studies [4-6]. 

At finite temperature, because quarkonium dissociation provides a probe of the pro- 
perties of the QCD medium created in heavy- ion collisions [7], Wilson loops and other 
observables related to Wilson lines have been widely used and measured on the lattice 
(see [8] for a review). The periodic boundary conditions that characterize the fields in 
a thermal field theory allow indeed for the definition of other gauge-invariant operators 
besides the Wilson loop. An example is the trace of the Polyakov loop, i.e. a Wilson 
line spanning the entire compactified time axis, and correlators thereof. Gauge-dependent 
quantities, such as the trace of the product of two Polyakov loops, have also been studied. 

In this paper, we focus on the cyclic Wilson loop, i.e. a rectangular Wilson loop whose 
time extent spans the entire Euclidean time axis, i.e. from to 1/T, where T is the tem- 
perature. This quantity can be seen as a gauge-invariant completion of the aforementioned 
product of two Polyakov loops. In particular, we study its renormalization. The moti- 
vation for this study comes from Ref. [9], where an order perturbative calculation in 
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dimensional regularization showed that the thermal expectation value of the cyclic Wilson 
loop is ultraviolet (UV) divergent after charge renormalization. 

It is known that a smooth Wilson loop is finite after charge renormalization in dimen- 
sional regularization [10, 11]. If the contour has cusps, additional UV divergences occur. 
These divergences are called cusp divergences and the coefficients multiplying them, and 
hence their anomalous dimensions, depend only on the angle at the cusp. Such divergences 
are renormalized through a multiplicative constant. This constant and the associated cusp 
anomalous dimension are known to two loops in QCD [12] and to three loops in = 4 
super symmetric Yang-Mills theory [13]. Furthermore, in [14], the intersection divergence, 
i.e. a UV divergence arising from an otherwise smooth contour intersecting itself, was first 
considered; it was shown that its renormalization is nontrivial and mixes all possible loops 
and correlators of loops sharing the same geometry. 

In this paper, we will show that the divergence appearing in the cyclic loop is an inter- 
section divergence, the intersection being caused by the periodic boundary conditions. The 
renormalization equation mixes the cyclic Wilson loop with the correlator of two Polyakov 
loops. This equation holds to all orders. We will explicitly check this up to order g^. 
Furthermore, the related renormalization group equations allow to resum logarithms asso- 
ciated with the intersection divergence; we will explicitly provide the resummed expression 
valid at leading logarithmic accuracy. 

The paper is organized as follows. In Sec. 2, we introduce some properties of Wilson 
loops and cyclic Wilson loops that will be relevant for the rest of the paper. In Sec. 3, we 
compute the cyclic loop at order g'^ at short distances r <C 1/T. Our result, UV divergence 
included, reproduces the short-distance limit of a calculation first presented in [9]. In Sec. 4, 
we introduce the renormalization of intersection divergences derived in [14] and apply it to 
the cyclic Wilson loop, obtaining its renormalization equation. The perturbative expression 
of the cyclic Wilson loop is used to determine the order g'^ contribution to the renormal- 
ization constant appearing in the renormalization equation. The renormalization-group 
equation is solved at leading logarithmic accuracy. In Sec. 5, we examine the structure of 
intersection divergences up to order g^; we show how our equation correctly cancels these 
divergences, providing a non-trivial test of the renormalization equation. In Sec. 6, we 
show that the same equation that renormalizes the cyclic Wilson loop at short distances 
renormalizes it at long distances. Our expression for the UV divergence at long distance 
disagrees with the corresponding one that can be found in [9]; we analyze the origin of 
the disagreement. Finally, in Sec. 7, we draw our conclusions, emphasizing the relevance 
of the result for a proper lattice evaluation of the cyclic Wilson loop. Some of the results 
presented here can be found in [15]. 

2 Definition and properties of the cycUc Wilson loop 

Before we turn to the issue of the renormalization of the cyclic Wilson loop, we will start 
by giving some general properties that will be instrumental in the following discussion. In 
Euclidean spacetime, the cyclic Wilson loop, Wc, is defined as the product of four straight 
Wilson lines, U, two of which extend in the time direction from to 1/T while the other 
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two stretch in the spatial direction r from — r/2 to 
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where (O) stands for the thermal average of the operator O. The Wilson lines, U , are 
defined as 
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where = A'^^T'^ is the gluon field in Euclidean spacetime with boundary condition 
A^{r,T + 1/T) = A^(r,r), V stands for path ordering of the colour matrices and Tr 
denotes the trace in colour space divided by the number of colours NcP' The four Wilson 
lines together give a rectangular path with sides 1/T and r. If we expand the cyclic Wilson 
loop in the coupling g and take heed of the path-ordering prescription, we get 
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The integral ffl dzj+i^ is performed along the rectangular path and ends at the point 

Zi as a consequence of the path-ordering prescription for the colour matrices. Note that 
the gauge fields are instead time ordered [3]. The thermal average (2.3) is a series of 
Feynman diagrams each of them characterized by n external points Zi lying on the contour 
of the rectangular path over which we integrate according to the path-ordering prescription. 
Sticking to the terminology found e.g. in Ref. [16], it is convenient to refer to these external 
points, each one of which gets a factor ig, as line vertices, as opposed to internal vertices 
connecting gluons to other gluons, ghosts or light quarks. The corresponding Feynman 
rule for line vertices reads 
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Figure 1 shows all order diagrams that can in principle contribute to the cyclic 
Wilson loop.'^ They represent the first non-trivial contribution in the perturbative series: 

VFf = [igfCp f dzip \j\\z2, D^, (zi, Z2) , (2.5) 

where Cp = Tr [T"T"] = (^N^ — Ij /{2Nc) is the quadratic Casimir in the fundamental 
representation. 

^Throughout the paper boldface characters wiU refer to three- vectors and italic characters to their 
modulus, e.g. r = |r|. 

^ We are working with the generators T" in = A'^T'^ in the fundamental representation of SU{Nc). 
^ This is the case for gauges in which Dot = 0, such as the Feynman or Coulomb gauges. In a generic 
gauge, there would be also four extra diagrams where the gluon connects temporal with spatial lines. 
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Figure 1. Diagrams tliat can in principle contribute at order as to the cyclic Wilson loop in gauges 
where Doi = 0. The time, r, runs along the horizontal axis and the spatial direction r along the 
vertical axis. 



It has been shown in Refs. [17] and [18] that the perturbative series for a Wilson loop 
can be rearranged and exponentiated as 



Wc = Y,C{-i)W{-i) = e^V 



762PI 



(2.6) 



where ^^(7) stands for the value of a diagram 7 without its colour factor given by C(7); 
(7(7) is the trace divided by Nc of all colour matrices, T", contracted with all colour 
structure constants appearing in the diagram 7. Equation (2.6) states that the sum over 
all diagrams can be exponentiated in such a way that the exponent contains the sum over a 
subset of these diagrams only, called two-particle irreducible (2PI) diagrams. However, the 
colour factors of the diagrams in the exponent have to be replaced by new factors C(7). 
These factors are called colour- connected [18] or maximally non-ahelian coefficients [17]. 
We will stick to the former expression and give an explicit definition below. 

The term two-particle irreducible was coined after the one-particle irreducible diagrams 
that occur e.g. in the resummation of the geometric series of self-energy diagrams for the full 
propagator (indeed, Ref. [16] uses the expression one-particle irreducible, but we will stick 
to the terminology found e.g. in Ref. [12]). Here "particle" refers to the contour only: if it 
can be cut in two places in such a way that the resulting pieces are not connected through 
gluons, then such a diagram is reducible. If there is no possible way to cut the contour 
twice that leads to two disconnected pieces, then such a diagram is called irreducible. This 
is shown in Fig. 2. The first diagram is reducible, because if cut as indicated, the two 
resulting pieces do not exchange gluons between themselves. These two pieces are given 
schematically on the right, dropping all specifics of the contour except for how the gluons 
are attached to it. The diagram shown in the second and third line cannot be separated 
like the first diagram, regardless of how it is cut. All possibilities reduce to one of the two 
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Figure 2. 

bottom). 



Example of a two-particle reducible (top) and an irreducible diagram (center and 



cases shown. Of course, if one would cut out a piece without any gluons attached to it, then 
every diagram would be reducible making this definition meaningless, so this trivial case 
is excluded. Sometimes 2PI diagrams are also called rainbow irreducible (cf. Ref. [12]). 

For further use it will be convenient to call a diagram connected, if every line vertex is 
connected to every other line vertex through gluons, internal vertices and possibly loops of 
quarks, gluons and ghosts, otherwise it will be called unconnected. Connected diagrams are 
irreducible, while unconnected diagrams, like those in Fig. 2, may or may not be reducible. 
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Figure 3. A list of all subdiagrams (called 7a etc.) obtained from the diagram on the left 70- 

With this terminology we can give a recursive definition of the colour-connected co- 
efficients. For a connected diagram, the colour-connected coefficient is equal to its colour 
factor. We can think of unconnected diagrams as combinations of connected parts. By 
combining two diagrams we mean putting their line vertices on the same contour. A subdi- 
agram of a given diagram can be obtained by removing any number of its connected parts. 
An example is given in Fig. 3: the four diagrams to the right represent all subdiagrams 
of the original diagram 70, obtained by removing one (in 7^ and 7^,) or two (in 7c and 
7(i) of the connected parts. In general there are several possible sets of subdiagrams that 
can be combined to form the original diagram. In our example these are {7a, 7d}, {76, 7c} 
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and {7c,7d,7d}- In order to obtain the colour-connected coefficient for an unconnected 
diagram we need to take the products of colour-connected coefficients of the subdiagrams 
for each of these sets and subtract them from the colour factor of the original diagram. If 
there are several ways of combining the subdiagrams to form the original diagram, then the 
product of colour-connected coefficients needs to be multiplied by the number of possible 
combinations. For example, there are two ways to combine subdiagrams 7^ and 7^: the 
left line vertex of 7^ can be put above or below the left line vertex of 7a; the same applies 
for the combination of {7c,7(i,7d}- If a subdiagram appears n times in a set, then the 
corresponding product needs to be divided by n\. In our example the colour-connected 
coefficient is therefore given by 

C{io) = C(7o) - 2 C(7a)C'(7d) - C(76)C(7c) - | C(7c) (^(7^))' • (2.7) 

The colour-connected coefficients of unconnected subdiagrams can be obtained by recursive 
application of this definition. The recursion ends when there are only connected subdia- 
grams involved. So e.g. for 7^ we have 

C(7a) = C(7a) - C(7c)C(7d) = C(7a) " C{^,)C{-fd) . (2.8) 

Although the definition does not require subdiagrams to be 2PI, we may neglect reducible 
subdiagrams in the calculation of colour-connected coefficients, since their colour-connected 
coefficients are zero. 




An illustration for Eq. (2.6) is given in Fig. 4. It shows a series of unconnected diagrams 
at O (og) with the colour factors written explicitly in front of each diagram. We have 



Cl = Tr 



rparparpbrj-ib 



and C2 = Tr 



rj-!arj-ibrj-iarj-it 



CfCa 



(2.9) 



where Ca = -^c is the quadratic Casimir of the adjoint representation. The ffist and last 
diagrams are reducible, the one in the middle is 2PI and its colour-connected coefficient 
is given by C2 = C2 — 2C|./2! = —CfCa/'^-, since is the colour-connected coefficient 
of a one-gluon diagram. The sum over the parts proportional to Cl of all three diagrams 
gives the first term in the second line of Fig. 4, which we interpret as the second-order 
expansion of the exponential of the 0{as) diagram, while the remaining 2PI diagram with 
its modified coefficient is interpreted as a term from the first-order expansion. 
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The exponentiation theorem greatly reduces the number of diagrams that we need to 
consider when calculating a Wilson loop, for we just need to consider 2PI diagrams.^ For 
the cyclic Wilson loop, we can neglect yet another group of diagrams. As Eq. (2.2) shows, 
Wilson lines are unitary operators, whose inverse is given by an otherwise identical Wilson 
line but with the direction of the contour integration reversed. Hence we can rewrite the 
cyclic Wilson loop as 



Wilson lines going from to 1/T on the time axis are Polyakov loop operators. Because 
they are related to the free energy of an infinitely heavy quark in the thermal medium [19], 
we will refer to them as quark lines. The other two Wilson lines will be called strings, 
because they act as gauge links between the two Polyakov loop operators. We will use the 
terms quark line and string to denote both the operators and the corresponding contours 
in spacetime. Because of the periodic boundary conditions on the Euclidean time, the two 
strings are at the same position and can be seen to combine into an adjoint Wilson line. If 
we expand one of the quark lines to zeroth order in then the strings cancel: 



In terms of diagrams, the zeroth-order expansion of a quark line corresponds to diagrams 
without any line vertices on that quark line. Equation (2.11) shows that in this case all 
diagrams with line vertices on the strings cancel against each other and what is left is a 
Polyakov loop. As the second line of the equation shows, this cancellation happens at any 
order in perturbation theory whenever the colour-singlet component of either of the two 
quark lines is considered. Some examples are given for illustrative purposes in Fig. 5. This 
leaves only diagrams with all line vertices on the same quark line or diagrams with some 
line vertices on both quark lines. Line vertices on the strings contribute only in the latter 
case. Because this cancellation is directly linked to the periodic boundary conditions and 
to the fact that the cyclic Wilson loop spans the entire compactified time axis, we will call 
it cyclicity cancellation. 

3 The cyclic Wilson loop at short distances up to O^a"^) 

Having illustrated some general properties of the cyclic Wilson loop, exponentiation and 
cyclicity cancellation, we can now turn to its calculation up to order g'^. In particular, we 
will investigate UV divergences. The analysis of this section will be extended up to order 
g^ in Sec. 5. 

*We remark that in an abelian theory all colour factors are one and only connected diagrams contribute 
to the exponent of Eq. (2.6). 



= (Tr [Uo (§, -§) UK_ (l/T, 0) C/J/^ (§, -§) Ur (l/T, 0)] ) 



(2.10) 



= (Tr [Uo (§, -§) uK (l/T, 0) U^^ (§, -|) U, (l/T, 0)] ) 
= [Uo (§, -§) (I + . . . ) U^o (§> -§) (l/T, 0)] ) 
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The cyclic Wilson loop is a thermal average, therefore, besides the scale 1/r that char- 
acterizes the correlation of two quark lines at a distance r, it depends on the temperature T. 
We will assume that these scales are much larger than the confinement scale, so that their 
contribution may be computed in perturbation theory. Other scales are also relevant: the 
in vacuum static energy, which in a weak-coupling regime is proportional to ag/r and the 
screening or Debye mass m^, which in a weak-coupling regime is proportional to gT. In 
order to be definite, we will assume l/r^T^m^)^ cts/r. The reason for this choice is 
that the same hierarchy was assumed for the correlator of two Polyakov loops in [20] and 
we will need the expression of the correlator of two Polyakov loops for renormalizing the 
cyclic Wilson loop at order g^. 

The cyclic Wilson loop is gauge invariant, so the choice of the gauge in which to perform 
the calculation is only a matter of convenience. A convenient gauge is the Coulomb gauge 
that we will adopt in this section if not stated otherwise, whereas the analysis of Sec. 5 
will not rely on a specific gauge choice. 

The 0{as) diagrams are shown in Fig. 1. The diagrams in the second row cancel 
because of cyclicity (see also Fig. 5). In the first row, only the diagram on the left, 
where the gluon connects to both quark lines, contributes to the cyclic Wilson loop. The 
other two diagrams vanish in dimensional regularization. The reason is that the contour 

rl/T 

integration along the time axis, i.e. / dr e*'^"'^ with cj„ = 2TTnT the bosonic Matsubara 

Jo 

frequencies, selects the zero mode n = 0, but then the remaining integration over the spatial 
momentum k depends neither on r nor on T and vanishes in dimensional regularization for 
being scaleless. Contributions from the scale rriD are of higher order. This line of argument 
is gauge independent, so the fact that the only contributing diagram is the first diagram 
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on the left in Fig. 1 is a gauge- independent statement. The diagram gives 

\nW^ = ^ + 0[al). (3.1) 
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Figure 6. All relevant diagrams at 0{al). As in Fig. 1 we restrict ourselves to gauges where 
Da^ = 0. 

At O (og) all irreducible diagrams that do not cancel through cychcity are shown in 
Fig. 6. The first three diagrams in the first row involve the gluon self-energy. This can be 
split up into a thermal and a vacuum part, where the thermal part is defined as the one 
depending on the Bose-Einstein or Fermi-Dirac distribution functions. In Coulomb gauge, 
the thermal part of the gluon self-energy at zero Matsubara frequency reads 



nS^(o,k) 
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(3.2) 



where nj is the number of massless fermions, and nsiq) = l/{exp{q/T) — 1) and np^q) = 
l/(exp(g/T) + 1) are the Bose-Einstein and Fermi-Dirac distributions respectively. We 
have taken the zero Matsubara frequency because, as before, this is the only contribution 
that survives the time integration. The gauge contribution to (3.2) can be read from [21], 
whereas the fermionic contribution can be found in textbooks such as [22]. The expression 
of IIqq (0,k) clearly shows that the thermal part of the self-energy is not UV divergent. 
Instead, all three diagrams that involve the gluon self-energy are IR divergent, but the IR 
divergences cancel each other [23] . Regarding the vacuum part of the gluon-self energy, its 
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expression in Coulomb gauge is known and can be read, for instance, from [24]. Summing 
together vacuum and thermal part, we get from the first three diagrams of Fig. 6 
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where /3o = llC^/3 — 2nj/3 and the first integral has been regularized in d = 3 — 2e 

- — r-T „ gives VOie), but, since the Oie) term multiplies 

(27r)" Airr 

a 1/e pole, we leave it uncomputed for the time being. The first line is the vacuum part 
and reproduces the well-known one-loop contribution to the QQ static potential [25, 26]. 
The UV divergence in the vacuum part can be removed by charge renormalization. In the 
second line, the term in square brackets simplifies to — nQT^(0, k)/k^ when k is integrated 
over the momentum regions k ~ 1/r and k T, whereas the part of IIqq (0,k) of order 
has to be kept unexpanded when integrating over the momentum region k ^ nif). 

The other two diagrams in the first row of Fig. 6 are finite; each contributes CFCA0i1/2. 
The first four diagrams in the second row give zero, because in gauges where time and space 
components do not mix, the three-gluon vertex for three longitudinal gluons vanishes. Also 
all five diagrams in the third row vanish in Coulomb gauge: the first four, because they 
involve scaleless integrals, and the last one because it is proportional to a function with 
vanishing support along the time axis. 

The last diagram in the second row, however, does give a contribution, which is ultra- 
violet divergent even after charge renormalization. The divergence comes from the vacuum 
part of the gluon connecting the strings. The thermal part is finite and will be given as 
a series expansion in rT, which, according to our adopted hierarchy of energy scales, is a 
small parameter. The radius of convergence for this expansion is rT < 1. We obtain 



(3.4) 

The previous calculation shows the specific advantages of the Coulomb gauge: UV 
divergences related to charge renormalization occur only in diagrams with gluon self-energy 
insertions, the uncancelled extra divergence arises only from the last diagram in the second 
row of Fig. 6, and several diagrams vanish. For comparison, in Feynman gauge only the 
diagrams with a three-gluon vertex in the second row vanish, while the two diagrams on 
the left of the third row are also divergent and contribute to both the divergence that is 
removed by charge renormalization and the one which remains after that. Since the cyclic 
Wilson loop is gauge invariant, the complete expression up to 0{a'^) is the same for both 
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gauges (we have explicitly checked this) and reads 
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(T) 

with IIqq (0,k) given by Eq. (3.2). Equation (3.5) agrees with the short-distance limit of 
the result that can be found in Ref. [9]; it thereby confirms their finding that the cyclic 
Wilson loop is not finite after charge renormalization. 

The divergence in the second line of (3.5) can be removed by charge renormalization. 
The divergence coming from the last diagram in the second row of Fig. 6, whose analytical 
expression is in the 1/e pole in the last line of (3.5), is very peculiar in that it is not of 
the form of a cusp divergence typically associated to a non-smooth Wilson loop. We recall 
that cusp divergences only depend on the angle 7 at the cusp. In Euclidean spacetime, the 
leading-order divergence reads in dimensional regularization^ [12] 



[1 -|- (vr — 7) cot 7] . 



(3.6) 



For four right angles, like in the case of a rectangular Wilson loop, Eq. (3.6) gives 
2asCi?/(7re); this divergence, being independent on r and T, can be removed by a multi- 
plicative factor. This is not the case for the divergence we are discussing here, which is of 
order and depends on the distance and on the temperature. In the next section, we will 
show how to properly renormalize this divergence. 



4 Renormalization 



In this section we will assume that charge renormalization has been carried out, so that 
we only need to concern ourselves with the remaining UV divergence in Eq. (3.5). This 
divergence is related to the specifics of the contour. Due to the periodic boundary con- 
ditions, the contour of a cyclic Wilson loop has support on a Euclidean cylindric space 
with d space dimensions extending from —00 to +00 and one compactified time dimension. 
Because the cyclic Wilson loop wraps around the full time dimension, the two strings run 
actually along the same line. This implies that the cyclic Wilson loop has no cusps, but 
only intersections, as shown in Fig. 7. 



^The angular dependence in a cut-ofF regularization was derived in [10]. 







Figure 7. The picture shows the contours of a non-cyclic (left) and a cyclic Wilson loop (middle). 
One can see how the cusp points turn into intersection points. The contour of the correlator of two 
Polyakov loops is shown on the right. 

It has been shown in Ref. [14] that the expectation values of Wilson loops with inter- 
sections cannot be renormalized by a single multiplicative constant. One has to consider 
instead sets of associated loops and loop correlators that mix under renormalization. By 
loop correlators we mean the expectation values (vacuum or thermal) of products of individ- 
ually traced loops. These sets of loops and correlators consist of all possible path-ordering 
prescriptions for contours that occupy the same points in spacetime and retain the same 
direction everywhere except at the intersection points. 

As an illustration consider the simple case shown in Fig. 8, which consists of a smooth 
curve intersecting itself once at a single point. Following the curve up to the intersection 
point, one can either go straight ahead thus following the rest of the contour (left picture in 
Fig. 8), or make a turn onto the way one has come from, thus splitting the contour into two 
separate loops (right picture in Fig. 8). Each of the two loops in the right picture, taken 
on its own, would have a cusp and be renormalizable through a multiplicative constant. 
However, when considering the product of the loops, there are new divergences coming 
from diagrams with gluon exchanges between the two loops. These new divergences are 
renormalized together with the smooth loop in the left picture, for which similar divergences 
arise at the intersection. More precisely, there exist linear combinations of the two loops 
shown in Fig. 8 that are finite and involve coefficients depending only on the angle at the 
intersection point. 



Figure 8. The two possible path orderings at the intersection point for a loop with one intersection 
point. In the figure on the right, the two loops are understood to be touching at the cusp points; 
the separation has been introduced to make clearer that one is dealing with two separate contours. 

In general, a Wilson loop may cross an intersection point several times and the angles 
at which the different lines enter the intersection point may all be different. In that case, 
the set of all associated loops is renormalized by a matrix of renormalization constants, 
which depend only upon the angles at the intersection point. When a loop has more than 
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one intersection point, then the set of associated loops takes on a tensor-hke structure with 
a renormahzation matrix for each intersection point. If there are additional cusps present, 
then those can be taken care of by multiplicative constants. So the general formula for the 
renormalized loops is^ 



(R) 



)Zi 



2J2 \ 



(4.1) 



Here the indices ik and jk label the different possible path-ordering prescriptions at the r 
intersection points, 9^ denote the sets of angles at those intersection points and (pi stand 
for the cusp angles at the s additional points of non differentiability. The loop functions, 
^jij2-jri '^^11 as the renormalized ones, W^^j^ v' ^'■^ defined such that there is a colour 
trace over each closed Wilson line, and each trace is normalized by the number of colours. 
The trace over closed Wilson loops ensures that all loop functions are gauge invariant. The 
coupling in W^^^^ is the renormalized coupling. The matrices Z are the renormahzation 
matrices. They are one at leading order in perturbation theory, while the explicit expression 
at higher orders depends on the adopted subtraction scheme. We will adopt the MS scheme. 

Now we want to apply the results of Ref. [14], which we have summarized above, and 
specifically Eq. (4.1) to the case of the cyclic Wilson loop. Although it may seem that the 
cyclic Wilson loop has a continuously infinite number of intersection points, namely all the 
points on the overlapping strings (see the second picture in Fig. 7), we need to care only 
about the two endpoints, for the Wilson loop contour does not lead to divergences in the 
other ones. As a consequence, we have four possible path orderings (two for each endpoint) 
when we consider the possible prescriptions at the intersection points, see Fig. 9. We will 
label the corresponding loop functions Wij, where i and j can assume the values or 1. 
These are explicitly given by 

:{l/T,0)U^, 



Woo 
Woi 

Wio 

Wn 



Tr 



I 
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Tr 
Tr 



Uo 
Uo 

ULil/T,0) 



^2' 



)Ur{l/T,0) 



Wr. 



v2' 



Tr 



2' 2 



)C/|(1/T,0) 
(1/T,0) 



Tr 



Pr 



UK{1/T,0) 



Uo (§, -§) UK{1/T,0)U^, (§, -§)] Tr \u^{l/T,0) 



Pr 



UK{1/T,0) 



Tr 



Uo 



2' ■ 



) Tr f/|(l/r,0) ) = Pc. (4.2 



We see that the four options give rise to two independent loop functions: the cyclic Wilson 
loop, Wc = Woo, and the correlator of two Polyakov loops separated by a distance r, 
Pc = Woi = Wio = Wii. We will call Pc simply the Polyakov loop correlator; its contour 
is shown in the rightmost picture of Fig. 7. 

We can represent Wij by a four component vector (VFoOi WqI) W^iOi W^ii) that gets 
renormalized by the tensor product of the two renormahzation matrices corresponding 



®The generalization to loop functions with more intersections performed in [14] relies on the assumption 
that the divergence structure at an intersection point is completely determined by the local characteristics 
of the contour at this point. This is certainly the case when for each intersection point there are at most 
two Wilson lines connecting it to other intersection points. For the cyclic Wilson loop we may therefore 
use Eq. (4.1), although it may not be applicable for more general loop functions. 
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Figure 9. The four different path orderings of the contour of the cyclic Wilson loop that correspond 
from left to right to the loop functions Woo, VVqi, Wiq and Wu respectively. The strings arc 
represented by the middle lines. 



to the two endpoints. Since the angles at both ends of the string are equal, also the 
renormalization matrices are equal. The renormalization equation reads 
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(4.3) 



Since the Polyakov loop correlator is finite, having neither cusps nor intersections, it holds 



that Pc 
on Wc- 



(R) 



( R) 

Pc- From this it follows that Ziq has to be zero, otherwise Pc would depend 



(4.4) 



The three equations involving Pc furthermore require Zu = 1 and Zqi = 1 — Zqo, which 
leaves only one independent renormalization constant Z = Zqq. So the renormalization 
condition for the cyclic Wilson loop and the Polyakov loop correlator reads 
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(4.5) 



We will now determine ZJ For the purpose of determining Z at 0{as), we just need 
to know that Pc = I +0{g'^). However, for the forthcoming analysis of Sec. 5, we will need 
Pc at 0{a^), which we give here. The expectation value of the Polyakov loop correlator is 
equal to the square of a single Polyakov loop, Pl, plus diagrams involving gluon exchanges 
between both loops. The one-gluon exchange diagram vanishes, because it is proportional 
to (Tr[r'^])^ = 0, so the first contribution comes from the exchange of two gluons [19, 27]. 
The result at O (og) for the same hierarchy of energy scales that we are assuming here, 
i.e. 1/r ^ T ^ mo ^ Oig/r, can be read from [20] or extracted in the appropriate limit 
from [19, 27, 28]. The Polyakov loop Pl = 1 + 6Pl has been computed in [9, 20]. The full 



^ Note that the renormalization constant Z is gauge independent for it relates gauge independent quan- 
tities . 
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O (og) result reads 

Pr 



(4.6) 
(4.7) 



1 2 y 
where = 



Let us now apply the renormalization equation (4.5) to the unrenormalized result (3.5). 
It is convenient to expand Z in powers of Og (understood as the 2e-dimensional coupling 
of dimensional regularization): 



Z=l + Zia.fi-'^' + Z2 (as/u"2")^ + 0{al) , 



where fi is the scale of dimensional regularization. To the purpose of fixing Z at order as, 
i.e. finding Zi, it is sufficient to write Eq. (4.5) as 



W]^) = ZWc + (1 - Z)Pc 



1 + 



rT 



T 



-irk 



+ Zia^ + 



(4.9) 



where the coupling is now the renormalized strong-coupling constant in the MS scheme 
and the dots stand either for finite terms at order a1 or for terms of higher order. Because 
Wc is finite, Eq. (4.9) fixes Zi; in the MS scheme it reads 



Ca fl 



IT \e 



-)E + In 47r 



CaI 

vr £ ' 



(4.10) 



where, for further use, we have defined l/e = 1/ e — + ln47r. 

After having removed the divergence from (4.9) according to (4.10), we can perform 
the Fourier transform of l/k^ in three dimensions, which gives l/(47rr), and write the final 
expression for the renormalized cyclic Wilson loop in the MS scheme: 



rT 



4tt 
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vr 



1 + 27ij - 2 In 2 + In /iV^ + ^ 



2(-l)X(2n),^^,2n 



rxy 
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(27r)2 



k2 + n£)(o,k) k2 



(4.11) 



The above expression is UV finite, the divergences having been reabsorbed either by the 
renormalization of the strong-coupling constant or by the subtraction of the intersection 
divergences along Eq. (4.5). 
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4.1 The cyclic Wilson loop at short distances at leading logarithmic accuracy 

Equation (4.11) is accurate up to next-to-leading order (NLO). It contains, however, loga- 
rithms in the renormalization scale ^ that may be potentially large. These logarithms can 

( R) 

be resummed by solving the renormalization group equations for Wc , which follow from 
the renormalization equation (4.5): 



"27r 



/3o 



(4.12) 



The renormalized coupling and the loop functions depend on the renormalization scale fi; 
the factor 7, 

1 (i „ (^s , 2\ 



j^-fi—Z = 2CA- + 0{at), 

Z dfl TT 



(4.13) 



is the anomalous dimension of the operator Wc 
group equation is trivial and reads 



(«) 



Pr. The solution of the renormalization 



(W^^'^-Pc) ifi) = (W^^^ - Pc) {1/r 



«s(/^) 

as(l/r) 



-4Ca//3o 



where we have made explicit the normalization scale dependence of Wc — Pc- Equation 
(4.14) implies that InVFc in the MS scheme at the scale /i may be written at NLO and 
leading logarithmic (LL) accuracy (i.e. including all terms of the type a^/ (rT) x (as In ^r)") 
as 
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(4.15) 



5 Renormalization up to 0{a^) 

The renormalization of the cyclic Wilson loop at 0{a'^) fixes the renormalization coefficient 
Zi, and only uses the fact that the Polyakov-loop correlator is l+0{g'^). In order to provide 
a nontrivial check of the renormalization equation, we will consider now the renormalization 
of the cyclic Wilson loop at 0{ag). At that order the expression (4.6) of the Polyakov loop 
correlator really matters. We will not attempt to compute the full expression of the cyclic 
Wilson loop at 0{a^), but we will just focus on its divergent contributions and on how 
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they are renormalized. For this purpose, we will start with a short discussion about the 
different types of divergences and their origin. We will summarize here the more detailed 
analyses found in Refs. [11] and [14]. Then we will focus on intersection divergences. In 
this section, if not otherwise stated, we will not rely on a specific gauge choice. 

5.1 Divergences 

Ultraviolet divergences come in general from integration regions in position space where two 
or more vertices are contracted to one point. In the case of internal vertices, one gets the 
usual divergences removed through charge renormalization. But for loop functions such as 
the cyclic Wilson loop, one also gets divergences from the contraction of line vertices along 
the contour. By contractions we mean only those that happen without moving vertices 
on a quark line to a string line or viceversa, and without altering the contour ordering of 
the vertices. So, for instance, in a diagram with one vertex on a quark line and two on 
the opposite, only the latter two can be contracted at the same point, either singular or 
smooth. For a generic diagram one has a superficial degree of divergence 

a; = l-7Vex, (5.1) 

at a smooth point (where the contour is differentiable and non-intersecting), and 

U = -iVex , (5.2) 

at a singular point (cusp or intersection); A'ex stands for the number of propagators con- 
necting the contraction point to uncontracted vertices.* There are therefore three possible 
situations that may lead to divergences related to line vertices: 

(1) all vertices are contracted to a smooth point, which leads to a linear divergence; 

(2) the contraction of vertices to a smooth point leaves an external propagator connecting 
a contracted to an uncontracted vertex: this leads to a logarithmic divergence that 
we will call line vertex divergence] 

(3) all vertices are contracted to a singular point, which gives a logarithmically divergent 
contribution; these are either cusp or intersection divergences. 



* The gluon propagator satisfies the periodic boundary condition: D^viQ, x) = D^,y{l/T, x). In the time 
interval < r < 1/T the thermal part of D^v{t,'x.) does not contribute to UV divergences, whereas these 
follow straightforwardly from considering the vacuum part of (r, x), for instance, in a covariant gauge ^, 



r{D/2 



1 + C 



^M- + (1 - C) 



D 



- 1 



where D — 4 ~ 2e is the number of dimensions. The counting rules (5.1) and (5.2) apply to covariant 
gauges. They also apply to the Coulomb gauge, where, however, some diagrams that would be divergent 
in a covariant gauge vanish (e.g. the first two in the last row of Fig. 6). The counting rules do not apply 
in general to singular gauges, where diagrams may exhibit a higher superficial degree of divergence. Since 
we are dealing with gauge-invariant quantities, we are allowed to exclude the case of singular gauges from 
our considerations. 
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Figure 10. Examples of linear divergences. The divergences arise when the vertices inside the 
dashed box are contracted to one point. 

Linear divergences are proportional to the length of the contour and can be removed by 
a factor that can be interpreted as a mass term; dimensional regularization removes these 
power-like divergences automatically [10]. Examples of diagrams with linear divergences 
are given in Fig. 10. In the notation adopted here, which follows the one in Ref. [14], the 
dashed box stands for integration regions where all vertices inside the box are contracted 
to one point. If the box includes a singular point, then the vertices are contracted to that 
point, otherwise they can be contracted anywhere inside the box. 

Line vertex divergences can be removed by using renormalized fields and couplings [11]. 
Cusp divergences arise from diagrams and integration regions as those depicted in Fig. 11. 
The one-loop divergence has been given in Eq. (3.6) as a function of the cusp angle 7. 
From it the renormalization constant for a non-cyclic Wilson loop (i.e. a Wilson loop with 
a time extension smaller than l/T) with four right-angled cusps can be inferred to be in the 
MS-scheme Z = exp [—2CFo:s^~'^'^/{iTe)] . Cusp divergences are absent in a cyclic Wilson 
loop. 




Figure 11. Contributions to a cusp divergence at 0{a^). 

We turn now to the intersection divergences of the cyclic Wilson loop, which are our 
main point of interest. They only appear when all vertices of a diagram or subdiagram are 
contracted to an intersection point. In all cases where at least one vertex is on the string, 
if every vertex of the diagram can be contracted to the intersection, then the contribution 
of the diagram cancels because of cyclicity. If all vertices are on a quark line, then the dia- 
gram contributes equally to the Polyakov loop, which is finite after charge renormalization. 
This leads to the conclusion that a connected diagram cannot give rise to an intersection 
divergence, because either all vertices can be contracted to an intersection point, in which 
case either the divergence cancels because of cyclicity or because it contributes to the 
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Polyakov loop, or it has at least one uncontr acted vertex and therefore it is finite according 
to Eq. (5.2). These different possibilities correspond to the three diagrams in the first row 
of Fig. 12. The first diagram has line vertices only on one quark line, so it also contributes 
to the Polyakov loop, which is finite. The second diagram has vertices on a string and on 
one quark line and thus cancels through cyclicity. The third diagram involves both quark 
lines, which means that the line vertices cannot be contracted to the same intersection 
point and therefore this diagram is finite according to (5.2). This exhausts all possible 
types of connected diagrams. 




Figure 12. The diagrams in the top row are connected and do not contribute to intersection 
divergences: they either cancel through cyclicity or contribute to the Polyakov loop. Conversely, 
some of the unconnected diagrams in the bottom row may show intersection divergences. See the 
text for details. 

Intersection divergences in a cyclic Wilson loop come from unconnected diagrams.^ In 
particular they come from unconnected diagrams made of at least one subdiagram with 
vertices on both quark lines and a subdiagram that is divergent once all its vertices are 
contracted to an intersection point. This can be understood by looking at the Feynman 
diagrams in the second row of Fig. 12. The first diagram is part of the Polyakov loop, which 
is finite, the second diagram cancels because of cyclicity, but the third and fourth diagrams 
are divergent, because we can contract the line vertices of the respective one-gluon and 
three-gluon subdiagrams to an intersection point. The periodic boundary conditions are of 
crucial importance here, because in order to have no external lines, the intersection point 
must be approached from the left and the right, as depicted in Fig. 13. 

As we just argued, diagrams contributing to intersection divergences need to contain 
at least one subdiagram with vertices on both quark lines. We call such subdiagrams 
bases. To be more specific, for each intersection-divergent diagram we define its basis as 
the subdiagram which is of highest order in as and still finite. In other words, the basis 
of a diagram is the subdiagram which is obtained by removing all parts that contribute to 
the divergence, but no more than that. Conversely, we can obtain all intersection-divergent 

^ Equation (2.6) allows us to restrict to 2PI unconnected diagrams. 
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Figure 13. Intersection divergences of subdiagrams. If tlie line vertices are contracted as in 
the left (to a smooth point) and middle diagram (to a singular point), we get respectively a line 
vertex divergence or a finite contribution because of the external line. But the periodic boundary 
conditions also allow for a contraction to a singular point without external lines, as shown by the 
right diagram. 

diagrams by considering all possible bases and combining them with all subdiagrams that 
lead to an intersection divergence. This classification of the intersection-divergent diagrams 
according to their bases and intersection-divergent subdiagrams provides the framework for 
the following analysis. 

All bases relevant for the computation of intersection divergences at O (a^) are given 
by the eight diagrams shown in Fig. 14. Note that bases may be reducible, provided that 
adding intersection-divergent subdiagrams makes them 2PI. In order to get all diagrams 
contributing to intersection divergences at this order, we need to add subdiagrams of O (q^) 
to the first basis, and a single gluon exchange to the other bases, in such a way that their 
line vertices can be contracted to at least one of the intersection points. 

The obtained diagrams may be divided into classes. A class of diagrams is made of 
diagrams that have the same basis and that can be transformed one into the other by just 
moving some of the line vertices across the intersection point from string to quark line or 
from quark line to string without changing their ordering along the contour of the Wilson 
loop. A class of O [a^) diagrams is shown in Fig. 15; a class will be typically represented 
by just one of its diagrams. 



5.2 Renormalization 

We consider first the effect of a single gluon when added to one of the diagrams in Fig. 14. 
It generates two classes of diagrams that lead to intersection divergences: one is shown in 
Fig. 16 with an intersection divergence at r/2, the other one would have the gluon attached 
to the lower quark like and an intersection divergence at —r/2. In Fig. 16, the contour 
integration for the left vertex of the gluon goes from (0, cr) to (a/T, r/2), and the contour 
integration for the right vertex goes from (—(1 — b)/T, r/2) to (0,dr), where we have used 
the periodic boundary conditions to shift the right vertex by —1/T in time; the grey area 
stands for a generic basis. Colour factors are considered separately. In a generic covariant 
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Figure 14. Bases relevant for the intersection divergences of O {c?^ 
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Figure 15. A class of 2PI diagrams represented by the diagram on the right. In Coulomb gauge 
the diagram on the right is the only one that does not vanish, see Sec. 3. 



jauge, the result reads: 
1 



vr 



^ + 27j;-ln4 + 2 + ln 7 ^ , ^\ +ln ^^ . ' ^\ ^r// 
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+ 2(1-0 In 



1 v(|-^)^'^' + «' 



i -c)^r2r2 + (1-6)2 



(l-6 + a)|c-(i|rr 



(5.3) 



The divergent part does not depend on the contour parameters a, 6, c, d, hence, it can be 
factored out from the contour integration. Because the contour parameters only appear 
in logarithms, they do not introduce new divergences when the contour integration for the 
basis is performed, since the basis is by definition free from intersection divergences. The 
fact that the divergence factorizes can be intuitively understood by noting that it comes 
from integration regions where all vertices of the added subdiagram are at the intersection. 
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so it is unaffected by the contour integration for the basis (see also Ref . [14] ) . Note that the 
divergent part of (5.3), i.e. —as/ir x 1/e, does not depend on the gauge fixing parameter ^; 
furthermore, the same divergent part occurs also in Coulomb gauge, in which case, however, 
only the integration along the strings contributes. The second class of diagrams with the 
gluon attached to the lower quark line gives the same divergence as in (5.3). 



Figure 16. Class of diagrams with a one-gluon subdiagram. 

At order O (og ) the only class of divergent diagrams is shown in Fig. 15. The divergence 
of that class is C2 x 2 x (— ag/vr x l/e)x (basis, i.e. the one-gluon exchange diagram, Os/r); 
C2 is the colour-connected coefficient of the diagrams, the factor 2 comes from the two 
intersection divergences at r/2 and — r/2, and the divergence is the one calculated in (5.3). 
The cancellation of this divergence in the MS-scheme fixes Zi to the value calculated in 
(4.10). 

Both Wc and Pc are equal to 1 at zeroth order in as- At that order the renormalization 
condition is automatically fulfilled: Z + {\ — Z) = \. A similar cancellation happens at 
any order (including odd powers of g from contributions of the Debye mass scale) when 
the renormalization constant multiplies diagrams that occur identically in both the cyclic 
Wilson loop and the Polyakov loop correlator. These are contributions to the cyclic Wilson 
loop coming from the gluon self-energy diagrams in the top row of Fig. 6 (second and third 
diagram) and from the diagrams in the bottom row of Fig. 6 (third and fourth diagram), 
and the contribution to the Polyakov loop correlator coming from the Polyakov loop, which 
is — 1 ~ 2 5Pl in Eq. (4.6). An example of cancellation is shown in Fig. 17. 



— 




>-v — 


■ Zx 







Figure 17. Cancellation of contributions from Wc and Pc- The Polyakov-loop correlator is repre- 
sented without strings. 

Let us discuss now divergences coming from O (ag) unconnected diagrams in the cyclic 
Wilson loop. We start by considering classes of diagrams whose basis is the second diagram 
in the first row of Fig. 14. The cancellation of the divergences carried by these diagrams 
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in the renormalization equation is a simple extension of the cancehation that happens at 
O (og) between the intersection divergence carried by the last diagram of the second row 
of Fig. 6 and Ziag times the one-gluon exchange diagram. This is illustrated in Fig. 18, 
where the bubble stands for any self-energy insertion. Also in this case, odd powers of 
g coming from the Debye mass scale, such as a term arising at order in our adopted 
hierarchy of energy scales, cancel out. 





finite 



Figure 18. The divergences of the two classes of diagrams cancel for the same reason as the lower- 
order classes of diagrams without self-energy insertions do. The giuon connecting the strings does 
not interact with the bubble. 

The cancellation of divergences carried by classes of diagrams whose bases are from the 
second row of Fig. 14 (shown in Fig. 19), which we will call collectively jt, is similar to the 
one discussed in the previous paragraph. They are cancelled by ZiUg times the respective 
basis diagram from the O {a^) expansion of Wc- As shown by Eq. (5.3) and the following 
discussion, adding one gluon to a basis gives a divergent factor, which is 2 x (— as/(7re)). 
The colour factor of the bases is iCfCa/'^, the colour-connected coefficient of diagrams 
like those in Fig. 19 is —iCpCj^/A. With this it is easy to see that the divergence of a class 
of diagrams with basis 7t is: 

(^-'-CfCI^ W (7t) = -^lOs (^ICfCa^ W (7t) , (5.4) 

where W (jt) denotes the value of a diagram belonging to "Jt without its colour factor. 
We see that all these divergences cancel in ZWc at O [a^)- 




Figure 19. Classes of diagrams with an intersection divergence and a three-gluon diagram as basis. 

With the cancellation of divergences associated to classes of diagrams whose basis is 
either the third or fourth diagram in the first row of Fig. 14, we come to the first non- 
trivial check of the renormalization equation (4.5). We call the basis with two ladder 
gluons (i.e. the third diagram in the first row of Fig. 14), -jjj, and the one with two 
crossed gluons (i.e. the fourth diagram in the first row of Fig. 14), ^x- Without colour 
factors the sum of these two bases gives just half of the square of the 0{as) diagram, 
so VF(7//) -|- W {'yx) = Q;g/(2r^T^). The colour-connected coefficients for the diagrams 
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with 7// and 7x as bases (these are the diagrams (d) and (e) of Fig. 20) are CfC\/4: and 
CfC\/2 respectively. The divergences associated to these classes of diagrams are of the 
type of Fig. 16 and bring a factor — as/vr x 1/e for each intersection point. In order to cancel 
these divergences, we have to collect also divergences coming from Zia^ x Wc that are not 
of the type cancelled by the Polyakov loop and divergences coming from the expansion of 
the exponential on the right-hand side of Eq. (2.6) (recall that the exponent is made by 
2PI diagrams only). We will highlight only the divergent terms and use the symbol "D" 
with the meaning "includes the divergent term". We have (compare with Fig. 20): 

(a) the renormalization constant ZiUs multiplying the second order expansion term of 

the one-gluon exchange diagram: Zia^ x Wc D ( ) x ^ ^ 



■K£ J 2 \ rT 

(b) the renormalization constant ZiUg multiplying the basis 'jx (whose colour-connected 
coefhcient is Ca = -CfCa/2): Zia^ x Wc D x (^-^^^W{jx] 

(c) the one-gluon exchange diagram times the intersection-divergent diagrams at O (og) 
from the expansion of the exponential: Wc D ^ « . . ^ « F s 



rT J \ TTE J \ rT 
(d) + (e) the divergent contributions from diagrams with bases 7// and jx- 



WcZ^i-'-^ 



\cFC\W(ciu) + \cfC\W{^x) 



Summing up these contributions, we get 



ZWc D - \ i ' [Wi^ii) + 2Wi^x) - Wijx)] - ir2^^ + 



Cl-\C.C,\^^. (5.5) 



2 J 2r2r2 vre 

This is exactly cancelled by the remaining contributions coming from the Polyakov loop 
correlator (4.6), given as (/) -1- (g) in Fig. 20: 



(l-Z)P.D^x 
ne 



„2 



2 2r2T2 



(5.6) 



The cancellation provides a non-trivial verification of the renormalization equation (4.5). 
Note that the cancellation of the diagrams in Fig. 20 is gauge independent, since Zi, 
the one-gluon exchange diagram, the Polyakov-loop correlator as well as the combination 
W (7//) -|- W (7x) are gauge invariant. 

Finally, we observe that the second-order expansion term of the renormalization con- 
stant times 0{as) diagrams of the cyclic Wilson loop gives a divergent term like 

ZWc D Z2al X ^ , (5.7) 

where Z2 may contain single or double poles in 1/e. This term is necessary to cancel diver- 
gent contributions coming from unconnected diagrams that have the one-gluon exchange 
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Figure 20. The sum of all these contributions is finite. It includes all diagrams with bases made 
of two gluons exchanged between the quark lines, and related contributions from both the cyclic 
Wilson loop and the Polyakov loop correlator (diagrams (/) and (g)). 
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diagram, i.e. the first diagram of Fig. 14, as a basis and subdiagrams of 0(a2) added. The 
different diagram classes are shown in Fig. 21. The sum of all divergent terms of these 
diagrams (after charge renormalization) determines the value of Z2. 

Because the right-hand side of Eq. (5.7) is proportional to l/{rT), in order to complete 
our proof of renormalizability at O (og) we have to make sure that the sum of diagrams 
in Fig. 21 only contains divergent terms proportional to l/{rT). The second and the third 
classes of diagrams in the first row and the classes of diagrams in the second row of Fig. 21 
do indeed fulfill this criterion. The reason is that, in order to get an intersection divergence 
from these diagrams, we have to contract all gluons but the one connecting the two quark 
lines to the intersection points, otherwise there would be a gluon propagator connecting the 
intersection point with an uncontracted vertex and from Eq. (5.2) we know that such a con- 
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figuration is convergent. As a consequence, since the divergence at the contraction does not 
depend on r, the complete r dependence comes from the one-gluon exchange, which is pro- 
portional to l/(rT). The same reasoning applies to the first and fourth classes of diagrams 
in the first row if we contract all gluons but the one connecting the two quark lines to the 
intersection points. Such kind of intersection divergence is called overall divergence in [14], 
to distinguish it from a subdivergence, which is a divergence that occurs when only one of 
the gluons is contracted to the intersection points. Consider the first diagram in Fig. 21. 
Its colour-connected coefficient is just like diagram (d) of Fig. 20. A subdivergence 

may come either by contracting the upper gluon connecting the strings to the intersection 
point at r/2 and keeping the finite part of the remaining diagram, or by contracting the 
lower gluon to the intersection point at —r/2 and keeping again the finite part of the re- 
maining diagram. The divergence may be read from (5.3) and is — as/C^re); the finite part 
is the same for the two situations and we call it W {'~f^)^^^^^'^\ Hence, the subdivergence of 

the first diagram in Fig. 21 is \cfCI ( ) VF(7+)(^°'*"). Because H^(7+)(fi'^it<=) has a 

4 \ 7re J 

more complicate functional dependence than just l/(rT), this is a divergence that is not 
canceled by the right-hand side of (5.7). Consider now the fourth diagram in Fig. 21 when 
none of the gluons is contracted to an intersection point. The colour-connected coefficient 

is C2 = -CfCa/2. The divergence is (ZiOs) x (^-^CfCa^ Ty(7+)(^"^*°) . Note that the 

finite part is the same multiplying the subdivergence in the first diagram. Summing up 
the (sub) divergences of the two diagrams we obtain 

ZW, D \cfCI W^(7+)(«'^^^^) + (-^) {-\cfCa) W^(7+)(«"^*^) = , (5.8) 

where we have used the value of Zi derived in (4.10). This completes our proof of renor- 
malizability of the intersection divergences at ©(as)- 

6 The cyclic Wilson loop at large distances 

Equation (4.5) establishes how the cyclic Wilson loop renormalizes, i.e. by mixing with 
the Polyakov loop correlator. Although we have tested the equation in a specific hierarchy 
of energy scales, i.e. 1/r » T » mD ^ Og/r, its validity is not bound to this hierarchy 
for it follows from general arguments based ultimately only on the UV behaviour of QCD 
and on the geometry of the loop functions [14]. In particular, Eq. (4.5) should also hold 
for cyclic Wilson loops at large distances, i.e. for rm/j ~ 1. 

Let us consider in the long-distance case rmo ~ 1 the UV divergences of the cyclic 
Wilson loop at 0{a^). To make contact with Sec. 3 we work in Coulomb gauge. The diver- 
gent diagram is again the last diagram in the second row of Fig. 6. The divergence comes 
from the vacuum part of the transverse gluon connecting the two strings. As expected, 
this quantity is insensitive to the low energy dynamics: the thermal part of the trans- 
verse gluon gives a finite contribution, the temperature providing a UV cut-off through the 
Bose-Einstein distribution, moreover transverse gluons are not screened. The temporal 
gluon connecting the quark lines is instead screened by the Debye mass. This is the only 
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difference witli tlie calculation performed in Sec. 3. Because the temporal gluon factorizes, 
at long distances the leading divergence of the cyclic Wilson loop due to intersection may 
be easily inferred from Eq. (3.4). Adding to it the leading-order diagram, which consists 
of the first diagram in Fig. 1, also with a screened temporal gluon, we get 

T 47rr T Attt e " ' ' ' ' ' 

where exp{—mDr) / (Airr) is the Fourier transform in three dimensions of the screened tem- 
poral gluon propagator, Z)oo(0,k) = l/(k^ + and the dots stand either for finite 
terms or for terms of higher order. It is straightforward to see that the expression in (6.1) 
is renormalized by (4.5) with the same renormalization constant Z computed in (4.8) and 
(4.10). This confirms, at least at leading order, our expectation that Eq. (4.5) renormal- 
izes the cyclic Wilson loop at any distance. Finally, we note that our conclusion and in 
particular Eq. (6.1) disagree with the long-distance finding of [9].^'^ 



7 Conclusions 



We have investigated the UV behaviour of the cyclic Wilson loop. Our short-distance, 
order-fl'^ calculation, summarized in Eq. (3.5), confirms the finding of Ref. [9], namely that 
the cyclic Wilson loop is divergent after charge renormalization and that its divergence is 
not a standard cusp divergence. 

In Sec. 4, we have shown how the periodic boundary conditions influence the diver- 
gences and the renormalization properties of the cyclic Wilson loop. The contour can be 
seen as having intersection points along the two strings, but only the two endpoints are rel- 
evant for renormalization purposes. Applying the intersection-divergence renormalization 
technique of Ref. [14], we have obtained the renormalization equation for the cyclic Wilson 
loop: Wc = ZWc + (1 — Z)Pc. This equation represents the main result of the paper and 
shows how the cyclic Wilson loop, mixes under renormalization with the correlator of 
two Polyakov loops, Pc, their mixing being determined by the renormalization constant Z. 
The equation holds for all SU{Nc) gauge theories; in the case of an abelian gauge theory, 
the cyclic Wilson loop and the correlator of Polyakov loops coincide and are finite. 

The disagreement may be traced back to the contribution from the non-zero modes to the transverse 
gluon connecting the strings in the last diagram in the second row of Fig. 6. In our computation, the 
relevant integral reads 

+7rrr ( - - 7B + In 47r + In ^ 
\e i 

The first line agrees with Eq. (A. 19) of [9], but the second line, which comes from the non- vanishing 
contributions from the double pole at r • q = in the two half-planes where the two Fourier exponentials 
separately converge, is missing. It is precisely the 1/e singularity in the second line that cancels the 
UV divergence coming from the zero modes, while the singularity in the first line, combined with other 
contributions, leads eventually to the result (6.1). We thank the authors of [9] for communication on this 
point. 
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We have determined the renormahzation constant, Z, up to order ^ obtaining Z = 
1 — CaOs/^"^^/'?'" X (1/^ — + ln47r) in the MS scheme. In Sec. 5, we have verified that this 
expression of Z reabsorbs all divergences of the type a^/{rT)'^; this is a non-trivial check 
involving the Polyakov-loop correlator at order a^. From the renormahzation constant we 
could extract the intersection anomalous dimension at one loop and solve the corresponding 
renormahzation group equations (4.12). The result provides the cyclic Wilson loop at LL 
accuracy (4.15), which is the novel computational outcome of this work. 

Finally, we observe that the renormahzation condition (4.5) is equivalent to stating that 
the combination Wc — Pc is multiplicatively renormalizable.^^ The combination Wc — Pc is 
therefore an ideal quantity to be computed on the lattice, while clearly the cyclic Wilson 
loop is not (see discussion in [9]). It could provide a new, independent and gauge- invariant 
lattice observable for the study of the thermodynamical properties of two static sources 
in a thermal bath, relevant for quarkonium physics in a quark-gluon plasma, being at the 
same time well suited for comparisons with analytic studies like the one performed in this 
work. 
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